Macroscopic description of current induced switching due to spin-transfer 
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We develop a macroscopic description of the current-induced torque due to spin transfer in a layered 
system consisting of two ferromagnetic layers separated by a nonmagnetic layer. The description is 
based on i) the classical spin diffusion equations for the distribution functions used in the theory of 
CPP-GMR, ii) the relevant boundary conditions for the longitudinal and transverse components of 
the spin current in the situation of quasi-interfacial absorption of the transverse components in a 
magnetic layer. The torque is expressed as a function of the usual parameters derived from CPP- 
GMR experiments and two additional parameters involved in the transverse boundary conditions. 
Our model is used to describe qualitatively normal and inverse switching phenomena studied in 
recent experiments. We also present a structure for which we predict only states of steady precession 
above a certain critical current. We finally discuss the limits of a small angle between magnetic 
moments of the ferromagnetic layers and of vanishing imaginary part of the mixing conductance. 

PACS numbers: 75.60.Ch,75.70.Cn,75.70.Pa 
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I. INTRODUCTION 



The magnetic moment of a ferromagnetic body can 
be switched without applying an external magnetic field, 
but only by transfer of electron spins carried in a spin 
polarized current. The concept of spin transfer has been 
introduced by SlonczewskiA and appears also in several 
papers of Berger 2 , Magnetic switching by a spin po- 
larized current has been now confirmed in extensive se- 
ries of experiments^*^. Spin transfer is also an im- 
portant turning point in spintronics. In spintronic phe- 
nomena of the first generation, like giant magnetoresis- 
tance (GMR)- or tunneling magnetoresistance (TMR) 8 , 
the magnetic configuration of a nanostructure is detected 
by an electrical current. On the contrary, in spin trans- 
fer experiments a magnetic configuration is created by 
a current. This possibility of back and forth magnetic 
switching opens new fields for spintronics. 

Current-induced magnetic switching (CIMS) has been 
clearly demonstrated 4 by experiments on structures 
F1/N/F2 consisting of two ferromagnetic layers of dif- 
ferent thicknesses separated by a nonmagnetic layer N. 
Starting from a parallel configuration of the magnetiza- 
tions in Fl and F2, a current exceeding a certain critical 
value can reverse the magnetic moment of the thinner 
magnetic layer to set up an antiparallel configuration. 
In turn, a current in the opposite direction can switch 
back the structure to the parallel configuration. With an 
applied field, the spin transfer mechanism can also gen- 
erate a steady precession of the magnetization, detected 
by oscillations of the current in the microwave frequency 
ranged. 



In the concept introduced by Slonczewskii, as well as 
in most theoretical modelsi2*ii*i£, the current-induced 
torque acting on the magnetization of a magnetic layer is 
related to the spin polarization of the current and, more 
precisely, to the absorption of the transverse component 
of the spin current by the magnetic layer. From CPP- 
GMR experiments one knows that the spin polarization 
of the current is due to spin dependent reflections at in- 
terfaces and to spin dependent scattering within the mag- 
netic layers. For CIMS, similarly, recent experiments 6 
have shown that the switching currents can be modified 
and even reversed by doping the magnetic layers with 
impurities of selected spin dependent scattering cross- 
sections. 

Both CPP-GMR and CIMS depend on spin accumu- 
lation effects. This is well known for CPP-GMR. For 
CIMS, this has been shown by experiments in which the 
spin accumulation profile is manipulated by introducing 
spin-flip scattering at different places in the structure^. It 
turns out that both the GMR effect and the spin transfer 
torque can be enhanced by introducing spin-flip scatter- 
ing outside a F1/N/F2 trilayer (in the leads) or reduced 
by spin-flip scattering in the nonmagnetic layer N 5 . This 
calls for a unified theory of CPP-GMR and spin trans- 
fer torque, taking into account spin accumulation, spin 
relaxation, and both interface and bulk spin dependent 
scattering. This is actually the direction of most recent 
theoretical developmentsAisiSii 3 ". 

The model we present in this paper fits directly with 
the interpretation of CPP-GMR data in the model of 
Valet and Fert (VF)s». Most of the parameters of our 
description can be derived directly from the analysis of 
CPP-GMR experimental data™, that is interface and 
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bulk spin asymmetry coefficients, interface resistance, 
layer resistivities, and spin diffusion lengths. As we 
will see below, two additional parameters, namely the 
real and imaginary parts of the mixing conductanc e) 13 ! 16 , 
are also needed. They can be derived from quantum- 
mechanical calculations 17 of the transmission of spin cur- 
rents at the interface under consideration. By introduc- 
ing into our model calculated values of the mixing con- 
ductance in addition to the set of parameters derived 
from GMR experiments, we calculate the current-induced 
torque for different types of structures in order to un- 
derstand how its sign, amplitude and angular variation 
depend on the spin asymmetry coefficients and spin accu- 
mulation effects. The calculations of our model are based 
on macroscopic transport equations similar to those de- 
rived from the Boltzmann equation approach of the VF 
modeli* for the CPP-GMR of multilayers with collinear 
magnetizations. Stationary charge and spin current are 
described by classical diffusion equations. We assume 
that the absorption of the transverse component of the 
spin currents is quasi-interfacial, as this has been jus- 
tified by a quantum description of the transmission of 
transverse spin current into a ferromagnetic layer. This 
assumption allows to derive some effective boundary con- 
ditions for the spin accumulation and spin current drops 
at the interfaces^, and this way also to calculate the 
torque acting on a magnetic film for an arbitrary an- 
gle between magnetic moments of the two ferromagnetic 
films in a structure. Thus, in a certain sense this extends 
to an arbitrary angle a small-angle description— that has 
been used for the interpretation of recent results 6 . 

The paper is organized as follows. Macroscopic equa- 
tions describing currents and spin accumulation inside 
the films are derived in section 2. The boundary condi- 
tions and general formulae for the torque in a four-layer 
structure are presented in sections 3 and 4, respectively. 
Numerical results for the structure with two magnetic 
films are presented and discussed in section 5. The lim- 
iting case of real mixing conductance is considered in 
section 6. The limit of a small angle between magnetiza- 
tions is discussed in section 7, whereas final conclusions 
are in section 8. 



A. Magnetic films 



The diffusion equation for arbitrary spin quantization 
axis takes then the formic 



b 9 2l = - 

dx' 



1 



2 



(1) 



where D is the diffusion 2x2 matrix in the spin space, 
1 is the 2x2 unit matrix, and r s f is the spin-flip re- 
laxation time. As it has been already mentioned in the 
introduction, we assume that the internal exchange field 
inside ferromagnetic metals is strong enough so that the 
component of the distribution function perpendicular to 
the local magnetization vanishes. Thus, the distribution 
function is diagonal when the spin quantization axis is 
parallel to the local spin polarization of the ferromag- 
netic system. Equation (1) can be then written as 



D 



dx 2 



1 



2t, 



D 



dx 2 



2t. 



sf 



-(/i-/t), 



(2) 



(3) 



where /f and /| are the distribution functions for spin- 
majority and spin-minority electrons, respetively. 

The above system of two equations can be rewritten 

as 



9 2 (/ T -/|) = 1 
I 



dx 2 



2 (/T - fl), 
sf 



d 2 {h+h) _J 2 {f ] -f [ ) 



dx 2 



dx 2 



where 



1 _ 1 / 1 1 

l 2 sf 2 \l 2 I, 



(4) 



(5) 



(6) 



with I 2 = DijT s f and I 2 = _DjT s /, and rj defined as 



II. CURRENTS AND SPIN ACCUMULATION 
INSIDE MAGNETIC AND NONMAGNETIC 
FILMS 



We assume the electrical current in the multilayer is 
carried by free-like electrons of equal concentrations in 
all the layers and without any spin polarization at equi- 
librium. The distribution function / inside the films is a 
2x2 matrix in the spin space, and its spatial variation 
can be described by the diffusion equation. We assume 
the distribution functions are uniform in the plane of the 
films, and vary only along the axis x normal to the films. 
Let us consider first ferromagnetic layers. 



V = 



(7) 



Equations (4) and (5) can be rewritten in terms of 
the electro-chemical potentials /If (fii) for spin-majority 
(spin-minority) electrons as 



<9 2 (MT ~Mi) _ 1 



dx 2 -frdh-hl 



dx 2 



dx 2 



(8) 



(9) 
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The above equations are equivalent to the equations de- 
rived from the Boltzmann equation approach by Valet 
and Feii^. 

Solution of Eqs (8) and (9) gives 



and 



P\ = (l+r?)[Aexp(o;/Z s /) 



■ B exp(-x/l sf )} + Cx + G, 



1*1 = {V~ l)[Aexp(x/l af ) 



(10) 



It is convenient to write the spin current in the matrix 
form as 



(21) 



with jo = (it+il) being the total particle current density, 
and j z = (j| — J j.) being the total ^-component of the spin 
current. Thus, one finds 



and 



1 



p(Ef) 



fa = -C(Di+D i ) l 



(22) 



■ B exp(-x// a /)] + Cx + G, 



(11) 



where A, B, C and G are constants to be determined 
later from the appropriate boundary conditions. 
The electro-chemical potentials can be written as 



with 



and 



h = Moi + g°z 



Po = (mt +M;)/2 



(12) 



(13) 



. 9 =(/2 T -/2 i )/2. (14) 
Thus, the explicit forms for £Lq and g are 

= v[Aexp(x/l sf ) + Bexp(-x/l sf )] + Cx + G, (15) 

and 

g = Aexp(x/l s f) + B cxp(— x/l s f). (16) 

For an arbitrary quantization axis the particle and spin 
currents are given by the 2x2 matrix j in the spin space 



(17) 



where p(Ep) is the density of states at the Fermi level 
per spin (per unit volume and unit energy). When the 
quantization axis is parallel to the local spin polarization, 
one finds 

j] = -D^C-— [Aexp(x/l sf ) - Bexp(-x/l sf )] , 



P(E F ) 



(18) 



' fa = -D l C+^- [Aexp(x/l sf ) -Bexp(-x/l sf )} , 

Isf 

(19) 



P(E F ) 
where 



D = 2- 



' D 1 + D i 



(20) 



P(E F ) 



U = -C{Pi-D i ) 



2D 

- —[Aexp{x/l sf ) -Bexp{-x/l sf )]. (23) 

Isf 

The particle current jo is related to the charge current 
Io via Iq = ejo, where e is the electron charge (e < 0). 
Thus, positive charge current (flowing from left to right) 
corresponds to negative particle current (electrons flow 
from right to left). 



B. Nonmagnetic films 

Solution of the diffusion equation for the distribution 
functions inside nonmagnetic films leads to the following 
equation 



P = Mol + g ■ & 



(24) 



where & — (a Xl <J yi <7 z ) and in a general case all the three 
components of g are nonzero. The general solutions for 
jlo and g have the forms 



fio — Cx + G, 



(25) 



g = Aexp(x/Z s /) + Bexp(-z/Z s /). (26) 
The spin currents arc then given by 



with 



and 



3 = 2O0I + j • or), 



j = -2CD 



P(E F ) 



(27) 



(28) 



1 2D 

— — j = - [Aexp(z/Z s/ ) -Bexp(-x/l sf )}, (29) 

P(J^F) hf 

where now = = D. Of course, all the constants 
may be different in different layers. 
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C. Rotations of the quantization axis 

Distribution function and spin current in the magnetic 
films are written in the coordinate system with the axis 
z along the local spin polarization. In turn, the formula 
given above for the distribution function and spin cur- 
rent inside nonmagnetic films have general form valid in 
arbitrary coordinate system. It is convenient, however, 
to write them in the system whose axis z coincides with 
the local quantization axis in one of the adjacent ferro- 
magnetic films. Since the magnetic moments of the two 
ferromagnetic films are non-collinear, it will be necessary 
to transform the distribution function and spin current 
from one system to another. Thus, if the solution for 
electrochemical potentials in a given coordinate system 
has the form (24), then the solution in the coordinate 
system rotated by an angle ip about the axis x is still 
given by Eq.(24), but with g replaced with g' given by 

g' x = 9x, (30) 
g'y = g y cos ip + g z sin <p , (31) 

9z = ~9y sin <P + 9z cos (p . (32) 

Similar relations also hold when transforming spin cur- 
rent j from one coordinate system to the other one. 



III. BOUNDARY CONDITIONS AND TORQUE 

To determine the unknown constants that enter the 
general expressions for electric current and distribution 
functions inside all the magnetic and nonmagnetic parts 
of any layered structure, we need to specify boundary 
conditions, which have to be fulfilled by the distribution 
function and currents at each interface. Such boundary 
conditions were derived by Brataas et ap* within the 
phenomenological description, and here we will make use 
of them. 

Charge and spin currents across the normal-metal- 
ferromagnet interface (called in the following interfacial 
currents), calculated on the normal-metal side in the co- 
ordinate system with the axis z along the local quantiza- 
tion axis in the ferromagnet, can be written as^: 

e 2 JO = (G T + G;)Gtf - j$) + (G T - G;)(gf -g»), (33) 

e 2 j z = (G T -G i )(^-/l A, ) + (G T +G i )(.gf -.gf), (34) 

e 2 Jx = -2Re{G Ti }gf + 2Im{G Ti }gf, (35) 

e 2 j y = -2Re{G u } 5 f - 2Im{G u } 5 f , (36) 



where g N (g F ) is the spin accumulation on the N (F) side 
of the N/F interface, Gj and G| are the interfacial con- 
ductances in the spin-majority and spin-minority chan- 
nels, and G|| is the spin-mixing conductance of the in- 
terface, which comes into play only in non-collinear con- 
figurations. It is worth to point that the above boundary 
conditions are valid when there is no spin-flip scattering 
at the interface. 

The boundary conditions can be specified as follows: 
(i) particle current is continuous across all interfaces (in 
all layers and across all interfaces it is constant and 
equal to jo), (ii) the spin current component parallel to 
the magnetization of a ferromagnetic layer is continuous 
across the interface between magnetic and nonmagnetic 
layers, and (iii) normal components (perpendicular to the 
magnetization of a ferromagnetic film) of the spin cur- 
rent vanish in the magnetic layer and there is a jump of 
these components at the interface between magnetic and 
nonmagnetic films, described by Eqs (35) and (36). The 
above boundary conditions have to be fulfilled at all in- 
terfaces. The number of the corresponding equations is 
then equal to the number of unknown constants, which 
allows one to determine the spin accumulation and the 
charge and spin currents. 

Since the perpendicular component of the spin current 
is absorbed by the magnetic layers, the corresponding 
torque t per unit square, exerted on a ferromagnetic film, 
can be calculated as 

T = |ax£-j±«), (37) 

where and j±r are the normal (to the magnetiza- 
tion) components of the spin current at the left and right 
interfaces of the magnetic film, calculated on the normal 
metal side of these interfaces. In the simple case (as in 
Fig. 1), where there is no additional magnetic layer out- 
side the F1/N/F2 trilayer and no transverse spin current 
at the outer edges of the trilayer, r is simply given by 
t = —hj±fj/2 for Fl and r = ?ij_Lz,/2 for F2, where j±r 
and j±L have to be calculated in the nonmagnetic spacer 
layer (N) at the right interface of Fl and left interface of 
F2 (left and right interfaces of N), respectively. 

IV. TORQUE IN A SPIN- VALVE STRUCTURE 

The structure F1/N/F2 under consideration consists 
of two left (thick) and right (thin) magnetic films, sepa- 
rated by a nonmagnetic layer. The thick magnetic film 
is assumed to be semi-infinite, so it also plays a role of 
the left lead. The thin magnetic film is followed by the 
right nonmagnetic lead, also assumed to be semi-infinite. 
Thickness of the nonmagnetic spacer layer is do, whereas 
of the thin magnetic film is di- Such a structure is shown 
schematically in Fig.l. Both ferromagnetic films are mag- 
netized in their planes, and magnetization of the thin 
layer is rotated by an angle <p around the axis x (normal 
to the films) as shown in Fig.l. Axis z of the coordinate 
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FIG. 1: Schematic structure of the system studied in this 
paper. The system consists of thick (Fl) and thin (F2) fer- 
romagnetic films, separated by a nonmagnetic (N) layer. The 
thick magnetic film (similarly as the right nonmagnetic lead) 
is assumed to be semi-infinite, while the thin nonmagnetic 
and ferromagnetic films have thicknesses do and cfe, respec- 
tively. The arrows indicate orientation of the net spin of the 
magnetic films, with (p being the angle between the spins. 



system is along the net spin of the thick ferromagnetic 
film (opposite to the corresponding magnetization). In 
both ferromagnetic films the local quantization axes are 
along the local net spin, while as the global quantization 
axis we choose the local one in the thick ferromagnetic 
film. According to our definition, charge current Iq is 
positive when it flows along the axis x from left to right, 
i.e., from the thick towards thin magnetic films (electrons 
flow then from right to left). 

The in-plane component t\\ of the torque acting on the 
thin magnetic film can be written as 



tii = a s x (s x S), 



(38) 



where s and S are the unit vectors along the spin polar- 
ization of the thin and thick magnetic layers, respectively. 
The parameter a is a function of the charge current Jo 
(not indicated explicitly in Eq.(38)). Equation (38) can 
be rewritten in the form 



a sm<p, 



(39) 



where the torque t v is defined in such a way that positive 
(negative) torque tends to increase (decrease) the angle 
if (<fi G (0, 27t)) between spin moments of the films. 
The torque t v can be calculated from Eq.(37) as 

h .. h . . . 

T v = -~3 y \x=d a = -7j {jz sin if + jy cosy) \ x=da , (40) 



an thin magnetic films. Comparison of Eqs (39) and (40) 
gives 



a = — 



fy\x=d = -i: Uz+ jy COt <p) \ x=do . (41) 



2sinc/? 



The out-of-plane (normal) component t± of the torque 
may be generally written as 



Ti = tsxS, 



(42) 



where the parameter b depends on I (not indicated ex- 
plicitly). It can be calculated from the formula 



. J x \x=a 



, Jx \x=d 7 



(43) 



where j x (J' x ) is the x-component of the spin current in 
the nonmagnetic thin film taken at the interface between 
the two thin films and written in the coordinate system 
of the thick (thin) magnetic films (j x = j' x according 
to Eq.(30)). From Eqs (42) and (43) follows that the 
parameter b is equal 



b= — 



2s'm(f 



3x\x=da 



h 



2 sirup 



Jx la;— do • 



(44) 



By taking into account Eqs (35) and (36) one can re- 
late the torque directly to the spin accumulation in the 
nonmagnetic film taken at the interface with the thin 
magnetic layer. The in-plane (Eq.(40)) and out-of-planc 
(Eq.(43)) torque components can be then rewritten as 

T v = [ Re i G u}9y +Im{G n }s4] \x=d 
= [ Rc { G n}(9y cos ip + g z sintp) 



and 



+lm{G u }g x ] 



T x = ~z [Re{G n }g' x - Im{G Ti }.g;] | x=do 



(45) 



= ^[Re{G n }g x -lm{G n } 
x(g y cos(p + g z sm<p)]\ x = do . (46) 

Similarly, the constants a and b can be related to the 
spin accumulation via the formulae 



« = — — [ Re i G tl}9 y + Im{Gu}5x] \x=d 

c sin 



R c{G n }{g v cot p + g z ) + Im{G n }-? 



9x 



smip 



=4,47) 



and 



where j' y and j z , y are the components of spin current in 
the nonmagnetic thin film written in the local system 
of the thin and thick magnetic films, respectively, and 
calculated at the very interface between the nonmagnetic 



b = 



e 2 sin ip 



[-Re{G n }g' x + lm{G n }g' y ] \ x=da 



MG U }^- + lm{G n }(g y cot p + g z ) 



snup 



4-48) 
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Equations (40,43) and (41,44) are the final formula for 
the torque components and the parameters a and 6, ex- 
pressed in terms of the spin currents. Alternatively, Eqs 
(45) to (48) are the corresponding formula expressed in 
terms of the spin accumulation. For numerical calcula- 
tions one can use either the former equations or equiva- 
lently the latter ones. 



V. NUMERICAL RESULTS 

For numerical calculations it is convenient to define the 
bulk and interfacial spin asymmetry factors for ferromag- 
netic films according to the standard definitions**, 



p Ta) = 2p*(lT/?) 



and 



R 



HI) 



2iZ*(lT7). 



(49) 



(50) 



where p^ and pi are the bulk resitivities for spin-majority 
and spin-minority electrons, respectively; i?| and R± are 
the interface resistances per unit square for spin-majority 
and spin-minority electrons, whereas and 7 are the bulk 
and interfacial spin asymmetry coefficients. The formula 
(49) will also be used for nonmagnetic films (with = 0). 
The conductances G-\ and G± (see Eqs (33,34)) are then 
G] = l/R-\ and G[ = 1/R[. The mixing conductance 
Gjj is generally a complex parameter with the imaginary 
part being usually one order of magnitude smaller than 
the real part. 

The key bulk parameters which enter the description, 
i.e., mean free paths and diffusion constants can be ex- 
pressed in a free electron model by the parameters de- 
fined in Eq. (49) and the relevant Fermi energy Ep. In 
numerical calculations we assume the same Fermi energy 
for both magnetic and nonmagnetic layers. The diffusion 
parameters f-m) are then calculated from the formulae 
(assuming free electron like model for conduction elec- 
trons), 



D Hl) = o VF ^Kl)> 



(51) 



where vf — y/2Ep /m e is the Fermi velocity of electrons, 
and the mean free paths AjQ) are 



m e VF 



ne- 



PHD' 



(52) 



with m e denoting the electron mass and n = 
(l/6ir 2 )(2m e EF/h 2 ) 3 / 2 being the density of electrons per 
spin. Apart from this p(Ep) (see Eq.(17)) is given 
by p(E F ) = (l/4Tr 2 )(2m e /h 2 ) 3 / 2 E F /2 . For such a de- 
scription (based on free electron like model) one finds 
Aj/Aj = (1 — 0)/(l + 0), and the parameter 77 defined by 
Eq. (7) is determined by via the simple relation 
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FIG. 



acting on the 



V = -P. 



(53) 



2: Normalized in-plane torque t v — 
thin ferromagnetic film due to spin transfer, calculated as 
a function of the angle tp for parameters typical for Co/Cu 
system, as described in the text. The insets indicate the sim- 
ilarity between the angular dependence of the parameter a of 
the expression tii = «sx(sxS) and that of the spin accumu- 
lation amplitude g (note that g is normalized to |Jo| whereas 
a to Jo). The other parameters are: Ef = 7eV, do = 10 nm, 
d-2, = 10 nm. 



For nonmagnetic films we use the same definitions, but 
now AjQ) and are independent of the spin orienta- 

tion (the corresponding is equal to zero). Finally, the 
spin diffusion lengths will be assumed as independent pa- 
rameters and will be taken from giant magnetoresistancc 
experiments. 

The parameters for the thick ferromagnetic film can be 
generally different from those for the thin magnetic film. 
Similarly, parameters corresponding to the two nonmag- 
netic components of the structure can also be different. 
In the following, however, we assume that the nonmag- 
netic spacer layer (N) and the right nonmagnetic lead are 
made of the same material. 

The following four different situations have been re- 
cently studied experimentally 6 . 

(i) fa = fa > 0, 71 = 72 > 0, which corresponds to 
Fl and F2 of the same material with positive spin asym- 
metries for both bulk resistivities and interfacial resis- 
tances (this means spin-majority electrons are less scat- 
tered both inside the layers and at the interfaces). 

(ii) (3i > 0, 71 > 0, 02 < 0, 72 < 0, which corresponds 
to different materials for Fl and F2, with positive spin 
asymmetries for Fl and negative spin asymmetries for 
F2. 

(iii) Pi =02 < 0, 71 =72 < 0, which corresponds 
to the same material for Fl and F2, with negative spin 
asymmetries for both bulk resisitivities and interfacial 
resistances. 

(iv) Pi < 0, 71 < 0, P2 > 0, 72 > 0, which corresponds 
to different materials for Fl and F2, with negative spin 
asymmetries for Fl and positive spin asymmetries for F2. 

One of the systems within the category (i) is Co/Cu 
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structure, that has been extensively studied experimen- 
tally. For the bulk resistivities and the interface resis- 
tances we take the experimental values obtained from 
the GMR measurements-—. Accordingly, for the Co lay- 
ers we assume p\ = p\ — 5.1 /if2cm, /3i = {3% = 0.51, 
ig^r = = 60nm, whereas for the nonmagnetic Cu lay- 
ers we assume Pq — 0.5 /ificm, l^y = 10 3 nm. 

In turn, for the Co/Cu interfaces we assume R\ = 
R* = 0.52 • 10" 15 fhn 2 and 71 = 72 = 0.76. In prin- 
ciple, the corresponding mixing conductance G-^ could 
be derived from the angular dependence of the CPP- 
GMR. However, in practice there is a large uncertainty 
in this derivation and, according to the experimentalists 
who have performed this type of experimen1ji2i2£i, there is 
no reliable experimental information on G^ from GMR. 
Therefore, we assume the value calculated in a free- 
electron model corrected by certain factors taken from 
ab-initio calculations by StilesSi. For free electron gas 
and no reflection at the interface (we assumed the same 
Fermi energy in all layers), one finds the following rela- 
tion between the spin current j' y and the spin accumu- 
lation g' y components (written in the coordinate system 
of the thin magnetic film and taken in the nonmagnetic 
spacer at the interface with the sensing (F2) magnetic 
film); j' y = p(E F )v F g' y /2 = 2G^g' y /e 2 . Here, G^ is 
the Sharvin mixing conductance in the limit when there 
is no reflection at the interface, G^ = e^kp/Airh, with 
the Fermi wavevector corresponding to the Fermi en- 
ergy Ep . Reflection from the interface can be taken into 
account effectively via a phenomcnological parameter Q, 
writing Re{G-f|} = QG^. In the case of Co/Cu system 
this factor is roughly equal to 0.925 according to Stiles 21 . 
Thus, in the following numerical calculations we assume 
Re{G u } = 0.542 • 10~ 15 flm 2 . As for the imaginary 
part, Im{G|^}, we determine it assuming the same ratio 
Im{G||}/Re{G||} as that following from ah initio calcu- 
lations. Thus, we assume Im{G^} = 0.016 • 10 -15 fim 2 . 

Numerical results on the in-plane and out-of-plane 
components of the torque as well as on the corresponding 
parameters a and b are shown in Fig. 2 and Fig. 3, respec- 
tively. The torque and the corresponding parameters a 
and b arc normalized to hIo/\e\. Within the linear model 
assumed here the spin accumulation and spin currents 
are linear functions of the charge current, so the curves in 
Fig. 2 are the same for arbitrary magnitude of the charge 
current Iq. Note that the sign of torque changes when Iq 
is reversed. 

Figure 2 implies that a positive current (Io > 0) 
tends to destabilize the parallel configuration and can 
switch it to antiparallel one above a certain threshold 
value. On the other hand a negative current tends to 
destabilize the antiparallel configuration. This behav- 
ior can be defined as a normal current-induced magnetic 
switching 6 . The numerical calculation of the switching 
currents is not the subject of our study here. We only 
note that, by using standard expressions of the switching 
currents as a function of the torque, magnetization, ap- 



0.003 




FIG. 3: Normalized out-of-plane torque, calculated as a func- 
tion of the angle ip for the parameters typical for Co/Cu sys- 
tem. The inset shows the corresponding parameter 6. The 
other parameters are the same as in Fig. 2. 



plied field, anisotropy field, Gilbert coefficient and thick- 
ness of the thin magnetic layer, one finds switching cur- 
rents of the order of 10 7 A/cm 2 . What we want to em- 
phasize is the physical picture emerging from the plots 
of Fig. 2. The main point is the similar angular depen- 
dence of the spin accumulation amplitude g = |g| in 
the nonmagnetic spacer and of the coefficient a in the 
expression (39) for the torque. This results from the 
relation between the transverse spin current and spin 
accumulation in the boundary conditions involving the 
mixing conductance, Eq. (36-37). The interpretation of 
the angular dependence is straightforward. The spin ac- 
cumulation g is larger in the antiparallel configuration, 
when the spin direction predominantly transmitted by 
the thick layer is slowed down by the thin layer, i.e. 
g(cp = 7r) > g(ip = 0). However, g(<p) does not increase 
monotonously from g(ip = 0) to g(ip = 7r); it begins with 
a decrease before getting at g(tp — tt). This initial de- 
crease follows from the enhanced relaxation of the spin 
accumulation due to the efficient pumping of transverse 
spins as ip departs from zero. This transverse pumping 
decreases to zero when ip tends to tt and g goes up to its 
maximum value g(ip = tt). The prefactor a — t v / sirup 
follows the same type of variation, with simply a slightly 
smaller initial decrease that can be explained by argu- 
ments related to the orientation of g, s and S (this will 
appear more clearly in section VI). Finally t v varies as 
asmtp, starting from zero as a(cp = 0)ip and going back 
to zero at ip = tt as a(ip = tt)(tt — ip), that is with a 
steeper slope if a((p — tt) > a(ip — 0). The sort of a 
shoulder seen in Fig. 2 for ip slightly below 7r/2 is related 
to the minimum in g and a at about this angle. 

The perpendicular component of the torque, that is 
the component coming from the imaginary part of Gfj., 
is shown in Fig. 3, and is rather small, much smaller than 
the in-plane component (it would vanish for Im{G||} = 
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0). Therefore, it plays a negligible role in the switch- 
ing phenomenon, at least for Co/Cu(lll). Although for 
other interfaces Im{Gf|}/Re{Gf|} is not as small as for 
Co/Cu(lll), the imaginary part of the mixing conduc- 
tance leads generally to perpendicular component of the 
torque which is definitely smaller than the in-plane one. 
Therefore, in the following we will deal only with the 
in-plane torque. 

Let us consider now the remaining three cases de- 
scribed above; (ii), (iii) and (iv). For simplicity, we as- 
sume the same absolute values of the bulk and interface 
spin asymmetry coefficients (3 and 7, but their signs are 
adapted to each situation. The other parameters are the 
same for all the cases. The torque corresponding to the 
four situations is shown in Fig. 4. 

The solid curve in Fig. 4(a) corresponds to the case (i), 
i.e. /3i = P2 > 0, 71 = 72 > 0. This curve is the same 
as that shown in Fig. 2, so that we will not discuss it 
here more. Let us change now the sign of the spin asym- 
metry parameters of the thin magnetic layer (case (ii) 
with /?2 < 0,72 < 0). The corresponding torque, shown 
by the dashed line in Fig. 4(a), has the same sign as in 
the case (i), so that the switching is still normal (io > 
destabilizes the parallel state) . With opposite spin asym- 
metry coefficients in Fl and F2, the spin accumulation 
g is larger in the parallel (tp — 0) state than in the an- 
tiparallel (tp = tt) one, so that the torque now starts from 
zero at tp — with a slope that is steeper than the slope 
corresponding to the return point to zero at tp = tt. 

The solid line of Fig. 4(b) corresponds to the case (iii) 
with negative values of all the spin asymmetry coeffi- 
cients. Compared to Fig. 4(a), the sign of the torque is 
now reversed, which means an inversion of the switching 
currents, as observed in FeCr/Cr/FeCr structures^. As 
the spin accumulation g in systems with similar materials 
for Fl and F2 is larger in the antiparallcl configuration, 
the slope is higher at the point where the torque comes 
back to zero at ip = tt. 

Finally, for the dashed line of Fig. 4(b) corresponding 
to the case (iv) (Pi < 0,71 < 0,/3 2 > 0,72 > 0), the 
switching is also inverse. This corresponds to the case 
studied in Ref.|(| with Fl (F2) and N corresponding to 
NiCr (permalloy) and Cu, respectively. Now, the spin 
accumulation is larger in the parallel state, so the torque 
approaches zero at ip = and p — it, with the slope 
larger in the former case (at tp = 0). 

Basic characteristics of the switching phenomena in all 
the four cases are gathered in Table I, where the sign 
of the torque for positive current (io > 0) is given for 
tp in the range < tp < tt. The normal/inverse switch- 
ing phenomenon is correlated there with the sign of the 
corresponding GMR effect. 

From the results described above follows that it is the 
spin asymmetry of the thick magnetic film, which de- 
termines whether the switching effect is normal or in- 
verse. When this spin asymmetry is positive (negative), 
one finds a normal (inverse) switching phenomenon. It 
is also interesting to note that when the spin asymme- 
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FIG. 4: Torque as a function of tp for the four situations 
described in the text, and for indicated values of the spin 
asymmetry parameters. The other parameters are the same 
as in Fig. 2 



tries of both magnetic films have the same sign, the 
structure shows normal GMR effect, whereas when they 
are opposite, the corresponding GMR effect is inverse 
as shown in many CPP-GMR measurements 22 . Exper- 
imental examples of the four behaviours of Table I can 
be found in the F1/N/F2 trilayers of Ref[6] respectively 
for NiFe/Cu/NiFe (i), NiFe/Cu/NiCr (ii), FeCr/Cr/FeCr 
(iii) and NiCr/Cu/NiFe (iv). 

For the parameters used in numerical calculations de- 
scribed above the current-induced torque vanishes in 
collincar configurations and one of them (either parallel 
or antiparallel) is unstable. This leads to either normal 
or inverse switching phenomena. An interesting situa- 
tion can occur when the amplitudes of spin asymmetry 
in the thick and thin layers are different. In Fig. 5, we 
show an example of torque calculated for Pi =0.1 and 
71 = —0.1, Pi = 0.51 and 72 = 0.76. Let us first con- 
sider the case with io > 0. The torque is positive when 
tp increases from zero, then comes back to zero at tp — tp c 
and becomes negative between tp c and tt. This means 
that the torque tends to destabilize both the parallel and 
antiparallel states. Above some threshold value of the 
current for the instability of the parallel and antiparallel 
states, the only solution is a steady state precession (in 
the absence of anisotropy and demagnetizing field, there 
would be a stable equilibrium at an intermediate orienta- 
tion between and tt in the layer plane, but the general 
solution is a precession). 

In trilayers which have been studied up to now, steady 
precessions with generation of microwave oscillations 
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Situation 


Spin asymmetries of magnetic layers 


Sign of torque (0 < 99 < 7r) 


Sign of GMR 


Thick layer 


Thin layer 


Torque 


Type 


(0 


(At)>o 


(At)>o 


r > 


normal 


GMR > (normal) 


(") 


09,7) >0 


(A 7) <0 


r > 


normal 


GMR < (inverse) 


(iii) 


(A 7) <o 


(A 7) <o 


r < 


inverse 


GMR > (normal) 


(iv) 


(At)<o 


(At)>o 


r < 


inverse 


GMR < (inverse) 



TABLE I: Characteristics of the current induced switching in the four cases studied in this paper. Correlation between the 
sings of the torque and GMR is also given. 



"7l = Pi = 0-1 0.8 
0.6 



/o<0 
0.04 \ 




FIG. 5: Torque for the case when the spin asymmetry factors 
of the thick ferromagnetic film, /3i = 0.1 and 71 = —0.1, are 
significantly different from those of the thin one, P2 = 0.51 
and 72 = 0.76. The torque is now normalized to |To| and the 
two curves correspond to To > and To < 0, as indicated. 
The angle tp c corresponds to the point where the torque due 
to spin transfer vanishes. Inset shows variation of the angle 
ip c with the spin asymmetry factor 71 for several values of j3i 
indicated on the curves. The other parameters are as in Fig. 2. 



have been generally observed when a magnetic field is 
applied and in a given range of current density^. The pos- 
sibility of obtaining microwave oscillations at zero field 
would be of great interest for several devices. Xiao et 
a&* have predicted that microwave oscillations can be 
observed at zero field and in a given range of current 
density when, for asymmetric structures, the spin trans- 
fer and damping torques have markedly different angu- 
lar dependences, so that their sum has a wavy angular 
dependence. In our structure of Fig. 5, the existence of 
steady precessions at zero field has a different origin. It 
comes from the wavy angular dependence of the torque 
itself. This structure, with Iq > would be of interest 
for the generation of oscillations at zero field and at any 
value of the current above some threshold value. 

The behavior of Fig. 5 for Io < has a different interest. 
Now, above some threshold value of the current, the spin 
transfer torque stabilizes both the parallel and antiparal- 



lel states of the trilayer. In other words, it increases the 
damping of the system in both configurations. This is of 
interest for some devices, for example for stabilizing the 
configuration of read heads based on CPP-GMR against 
fluctuations or, at least, avoiding spin transfer induced 
fluctuations. 



VI. DISCUSSION IN THE LIMITING CASE OF 
REAL MIXING CONDUCTANCE 

As it has been already mentioned before, the imaginary 
part of the mixing conductance is usually small. When 
Im{G||} = 0, the formulae (45) and (46) for the torque 
components acquire a simpler form. First of all, the out- 
of-plane components of the spin accumulation and spin 
current vanish, g x = and j x — 0. The in-plane torque 
can be then written in the form 

U = 

= -^ G u(9yCOSf + g z sm(p)\ x=do , (54) 

where G-\[ is real. In turn, the out-of-plane torque van- 
ishes then exactly, 



7~x — 0* 



(55) 



When the spin accumulation in the layer N at its in- 
terface with the layer F2 forms an angle <p g with the axis 
z, then one finds g' = gsin(ip — tp g ) and the torque may 
be written in the form 



r G n gsin(</j - ip g ), 



(56) 



where g is the absolute value (amplitude) of the spin 
accumulation at the interface. According to Eq.(41) the 
parameter a may be then expressed in the form 



h sin(<p - (p g ) 

-^U5 : : 

e z sin (p 



(57) 



The formulae (56) and (57) relate the torque and the 
parameter a to the amplitude g of spin accumulation 
in the layer N at its interface with F2, and to the si- 
nus of the angle ((p — (p g ) between the spin accumula- 
tion g and the polarization axis of the thin ferromag- 
netic layer F2. It is interesting to look at the angu- 
lar variation of these parameters. We have calculated 
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FIG. 6: Amplitude of spin accumulation at the interface (a), 
the angle tp g (solid line in (b)), the factor sin(yj — tp B )/ simp 
(dashed line in (b)), the torque r v (dashed line in (c)), and 
the parameter a (solid line in (c)). All curves are calculated 
for ImGji = and for positive current, To > 0. The other 
parameters are as in Fig. 2. 



before the upturn to g in the antiparallel state is due to 
the relaxation enhancement generated by transverse spin 
pumping in a non-collinear state. As for the orientation 
of g, one can see from the variation of tp g in Fig. 6(b) 
that the vector — g has first an intermediate orientation 
between those of S and s, following the rotation of s at an 
angle of about 0.7<p from s. Then — g comes back to the 
orientation of S and, as tp tends to n, the angle between 
— g and S tends to as — 0.55(7r — tp). The result of this 
angular variation of the spin accumulation is that the fac- 
tor — sin(<p — <p g )/ sin tp does not change significantly and 
varies only between about 0.52 and 1.05 for the system 
we have considered, as shown in Fig. 6(c). This indicates 
that the variation of the torque prefactor a, shown in 
Fig. 6(d), is mainly controlled by the variation of the spin 
accumulation amplitude g, with a small additional in- 
fluence of the factor — sin(<^ — ip g )/ shop (this influence 
explains, for example, that the dip around tp = 0.457T is 
less pronounced for a than for g). Finally the variation 
of t v with the angle tp, dashed curve in Fig. 6(d), reflects 
the variation of the product of a times sin tp. 

We have seen in section V that similar angular vari- 
ation of g and a = r»/ sirup also occurs in the gen- 
eral case with nonzero but small imaginary part of Gfj.. 
The explanation is the same as above, but only a lit- 
tle more complex because g has also a normal com- 
ponent. The general conclusion is that the torque is 
closely related to the spin accumulation in the nonmag- 
netic layer N. We put forward two practical consequences 
of this correlation. (1) The ratio between the switch- 
ing current amplitudes of the parallel— >antiparallel and 
antiparallel— ^parallel transitions, which reflects approxi- 
mately the ratio between the derivatives dr^/dtp at tp = 
and (p = 7r, simply reflects the ratio between g in the par- 
allel and antiparallel configurations. In particular this 
ratio is inverted when the sign of the spin asymmetries 
is inverted in one of the magnetic layers, as in the situ- 
ations (iii) and (iv) in section V. (2) More generally, the 
torque amplitude can be enhanced or reduced by enhanc- 
ing or reducing the spin accumulation in the nonmagnetic 
spacer. This has been confirmed by the experiments of 
Rcf. 6] in which the torque could be enhanced or reduced 
by introducing spin-flip scatterers at different places in 
the structure. 



them for the case (i) and in Fig. 6 we show the angu- 
lar variation of the spin accumulation amplitude g, the 
angle <p g: the factor — sm(<p — tp g )/ simp of Eq.(57), the 
torque amplitude t v , and the prefactor a of the expres- 
sion tii = asx(sxS). All the curves are shown there for 
positive current, I > 0. The amplitude g of the spin ac- 
cumulation, see Fig. 6(a), goes from its small value in the 
parallel state to its higher value in the antiparallel state, 
as expected when Fl and F2 have the same spin asym- 
metries. As discussed in section V, the initial decrease 



VII. COMPARISON WITH THE PREVIOUS 
CALCULATION OF THE TORQUE FOR tp 
CLOSE TO OR i 

A simple expression of the current-induced torque has 
been derived by Fert et aii& in the small angle limit, that 
is when the angle tp between the two magnetic layers {9 
in the notation of Ref . |l8j ) is small or close to n. This 
expression involves also parameters derived from CPP- 
GMR data and has been recently fitted with experimen- 
tal results of CIMS on samples doped in several ways&. It 
is interesting to compare the torques at small angle cal- 
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culated with this expression and by using the model of 
the present article. For simplicity, we will compare the 
two calculations when the mixing conductance is real, 
which allows us a clearer physical picture of the differ- 
ences between the two approaches. With zero imaginary 
mixing conductance (e = in the notation of Ref . [Tsl] ) . 
the expression of the torque at small angle is given by 
Eq. (5) of Ref. [3, 



P(AP) 

v F m N y 



P(AP) ■ 

Jm,N 



(1 



-d /\\ 



Actually the factor 1/2 in Eq.(61) comes from the val- 
ues ip g = n + ip/2 for i p in finitcsimally small (Q m = 9/2 
in the notation of Ref.[l]|) and (p g = it — (it — ip)/2 for 
7r — ip infinitcsimally small. These values of ip g are de- 
rived from a transverse spin conservation condition with 
the assumption of a constant orientation of g in a thin 
enough nonmagnetic layer. Our calculations in this pa- 
per show that, for a Cu layer of 10 nm, this assumption 
does not strictly hold and that the factor 0.5 must be 
replaced by 0.7 and 0.52 for <p close respectively to and 



P(AP) "> 
VFm F \ .P(AP) 
A ' 3m, Fl 



,-do/X 



s x (s x S), (58) 



where m^j AP ^ {j^jf^ ) ^ s ^ ne s P m accumulation density 
(spin current ) in the N layer at the N/F2 interface 
in the parallel (P) and antiparallel (AP) configurations, 
w p(ap)^.p(ap)-j ^ e same quantities in Fl at the 
Fl/N interface (m and j m are defined as positive for po- 
larizations in the majority spin direction of Fl, note also 
that, in the notation of Ref. 18|, one electron counts for 
1/2 in m and j m ) and A is the mean free path in N. 
The torque we have calculated in the preceding sections 
corresponds to the first and dominant term of this ex- 
pression in the limit do > > A (we will come back later to 
the meaning of the other terms). Omitting the last three 
terms, Eq.(58) can be written as 



P(AP) 

v F m N y 



s x (s x S). 



(59) 



The spin accumulation density m can be written as a 
function of the spin accumulation g of our paper in the 
following way: 

m e kp 

m = ^~T^ 9 ( 60 ) 
2-K l n 

so that Eq.(59) can be written as 

t = ^Gftg s x (s x S), (61) 
to be compared with the torque found in section VI: 



h sxn(<p 

--o^tld 

e sm ip 



<Pg) 



s x (s x S). 



(62) 



The only differences between Eqs.(61) and (62) are: 

(i) The replacement of G^ by G^. In section V, we 

assumed Gfj. equal to 0.925 xG^j 1 for Co/Cu interface. 
The factor 0.925 is also the factor t which wa s ap proxi- 
mated by 1 between Eq.(4) and Eq.(5) in Ref-filf). 

(ii) The replacement of the factor 1/2 by the value of 
— sin((y3 — ifg)/ sin(yj) for tp close to and n - approxi- 
mately by 0.7 and 0.52 in the case illustrated in Fig. 6(b). 



We thus conclude that the torques expressed respec- 
tively by Eq.(61), that is derived from the first term in 
the expression of the torque in Ref.[l8|. and by Eq.(62) 
derived in this article, differ only by a numerical factor 
not very different from unity. This factor will tend to 
unity for thinner N and will depart more from unity for 
thicker layer N. 



It remains to discuss why our calculations do not in- 
clude the three last terms of Eq.(58). This comes from 
the mixing conductance approximation^, that is the ap- 
proximation of the boundary equations for the transverse 
components, Eqs. (35,36). These equations express the 
diffusion transverse spin current generated by the dis- 
continuity between the finite value of gj_ (g_i_ denotes 
the spin accumulation component normal to the magne- 
tization) in N and its zero value in F. The finite value of 
gj_ is taken just at the interface, which assumes that the 
gradient of gj_ and the resulting variation of gx on a dis- 
tance of the order of the mean free path can be neglected. 
The second term of Eq.(58) takes into account the con- 
tribution from this gradient to the diffusion current. In 
addition, if N is thinner than the mean free path, a cer- 
tain amount of the diffusion current comes directly from 
Fl, which gives rise to the last two terms of Eq.(58). In 
conclusion, the calculation of the torque at small angle 
in this article has the advantage of a more accurate de- 
termination of the orientation of the spin accumulation 
in the nonmagnetic spacer, that is, for example, a more 
accurate determination of the factor sin(</? — ip g )/smip 
involved in Eq.(62). In the system we considered, the dif- 
ference is relatively small (for a 10 nm thick spacer, 5% 
and 40% for <p close respectively to and n), and it should 
decrease (increase) for thinner (thicker layers). On the 
other hand, with the boundary conditions of Eqs. (35, 36), 
we are not able to calculate the contributions to the dif- 
fusion current due to the gradient of spin accumulation 
and to the direct diffusion from the thick magnetic lay- 
ers. These contributions, corresponding to the last three 
terms of Eq.(58), must be taken into account for thin 
nonmagnetic spacers - the diffusion from the thick layer 
vanishes only when the spacer is thicker than the mean 
free path. We will introduce them in a further extension 
of our model. 
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VIII. CONCLUDING REMARKS 

In conclusion, we have presented a model of CIMS 
which is partly based on the classical transport equations 
derived from the Boltzmann equation standard model of 
CPP-GMR 14 . Additional boundary equations based on 
the concept of mixing conductance are used to describe 
the interfacial absorption of transverse spin currents in 
non-collinear magnetic configurations. This model ap- 
plied to Co/Cu/Co trilayers allows us to calculate the 
spin transfer torques as a function of the usual param- 
eters derived from CPP-GMR measurements (interface 
spin asymmetry coefficients and resistances, bulk spin 
asymmetry coefficients and resistivities, spin diffusion 
lengths), and the mixing conductance coefficients derived 
from ab initio calculations. We have also shown that 
the torque and its angular dependence is closely related 
to the spin accumulation in the nonmagnetic spacer and 
to its angular dependence. Enhancing the spin accumu- 
lations seems to be the way to reduce the critical cur- 
rents. The model has been also applied to situations 
with different spin asymmetries in the two magnetic lay- 



ers of F1/N/F2 structures to reproduce the inversion of 
the switching current recently obtained by reversing the 
spin asymmetry of the thick magnetic layer. By applying 
it to asymmetric structures, we have shown that steady 
precessions can be obtained in zero field. We have also 
pointed out some application limits for boundary condi- 
tions and indicated that certain corrections can be antic- 
ipated by going beyond these limits. 
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